3.29 \(\int \frac{c+d x}{(a+a \coth (e+f x))^3} \, dx\)

Optimal. Leaf size=183 \[ -\frac{c+d x}{8 f \left (a^3 \coth (e+f x)+a^3\right )}+\frac{x (c+d x)}{8 a^3}-\frac{11 d}{96 f^2 \left (a^3 \coth (e+f x)+a^3\right )}+\frac{11 d x}{96 a^3 f}-\frac{d x^2}{16 a^3}-\frac{c+d x}{8 a f (a \coth (e+f x)+a)^2}-\frac{c+d x}{6 f (a \coth (e+f x)+a)^3}-\frac{5 d}{96 a f^2 (a \coth (e+f x)+a)^2}-\frac{d}{36 f^2 (a \coth (e+f x)+a)^3} \]

[Out]

(11*d*x)/(96*a^3*f) - (d*x^2)/(16*a^3) + (x*(c + d*x))/(8*a^3) - d/(36*f^2*(a + a*Coth[e + f*x])^3) - (c + d*x
)/(6*f*(a + a*Coth[e + f*x])^3) - (5*d)/(96*a*f^2*(a + a*Coth[e + f*x])^2) - (c + d*x)/(8*a*f*(a + a*Coth[e +
f*x])^2) - (11*d)/(96*f^2*(a^3 + a^3*Coth[e + f*x])) - (c + d*x)/(8*f*(a^3 + a^3*Coth[e + f*x]))

________________________________________________________________________________________

Rubi [A]  time = 0.205172, antiderivative size = 183, normalized size of antiderivative = 1., number of steps used = 11, number of rules used = 3, integrand size = 18, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.167, Rules used = {3479, 8, 3730} \[ -\frac{c+d x}{8 f \left (a^3 \coth (e+f x)+a^3\right )}+\frac{x (c+d x)}{8 a^3}-\frac{11 d}{96 f^2 \left (a^3 \coth (e+f x)+a^3\right )}+\frac{11 d x}{96 a^3 f}-\frac{d x^2}{16 a^3}-\frac{c+d x}{8 a f (a \coth (e+f x)+a)^2}-\frac{c+d x}{6 f (a \coth (e+f x)+a)^3}-\frac{5 d}{96 a f^2 (a \coth (e+f x)+a)^2}-\frac{d}{36 f^2 (a \coth (e+f x)+a)^3} \]

Antiderivative was successfully verified.

[In]

Int[(c + d*x)/(a + a*Coth[e + f*x])^3,x]

[Out]

(11*d*x)/(96*a^3*f) - (d*x^2)/(16*a^3) + (x*(c + d*x))/(8*a^3) - d/(36*f^2*(a + a*Coth[e + f*x])^3) - (c + d*x
)/(6*f*(a + a*Coth[e + f*x])^3) - (5*d)/(96*a*f^2*(a + a*Coth[e + f*x])^2) - (c + d*x)/(8*a*f*(a + a*Coth[e +
f*x])^2) - (11*d)/(96*f^2*(a^3 + a^3*Coth[e + f*x])) - (c + d*x)/(8*f*(a^3 + a^3*Coth[e + f*x]))

Rule 3479

Int[((a_) + (b_.)*tan[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[(a*(a + b*Tan[c + d*x])^n)/(2*b*d*n), x] +
Dist[1/(2*a), Int[(a + b*Tan[c + d*x])^(n + 1), x], x] /; FreeQ[{a, b, c, d}, x] && EqQ[a^2 + b^2, 0] && LtQ[n
, 0]

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 3730

Int[((c_.) + (d_.)*(x_))^(m_.)*((a_) + (b_.)*tan[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> With[{u = IntHide[(a
+ b*Tan[e + f*x])^n, x]}, Dist[(c + d*x)^m, u, x] - Dist[d*m, Int[Dist[(c + d*x)^(m - 1), u, x], x], x]] /; Fr
eeQ[{a, b, c, d, e, f}, x] && EqQ[a^2 + b^2, 0] && ILtQ[n, -1] && GtQ[m, 0]

Rubi steps

\begin{align*} \int \frac{c+d x}{(a+a \coth (e+f x))^3} \, dx &=\frac{x (c+d x)}{8 a^3}-\frac{c+d x}{6 f (a+a \coth (e+f x))^3}-\frac{c+d x}{8 a f (a+a \coth (e+f x))^2}-\frac{c+d x}{8 f \left (a^3+a^3 \coth (e+f x)\right )}-d \int \left (\frac{x}{8 a^3}-\frac{1}{6 f (a+a \coth (e+f x))^3}-\frac{1}{8 a f (a+a \coth (e+f x))^2}-\frac{1}{8 f \left (a^3+a^3 \coth (e+f x)\right )}\right ) \, dx\\ &=-\frac{d x^2}{16 a^3}+\frac{x (c+d x)}{8 a^3}-\frac{c+d x}{6 f (a+a \coth (e+f x))^3}-\frac{c+d x}{8 a f (a+a \coth (e+f x))^2}-\frac{c+d x}{8 f \left (a^3+a^3 \coth (e+f x)\right )}+\frac{d \int \frac{1}{a^3+a^3 \coth (e+f x)} \, dx}{8 f}+\frac{d \int \frac{1}{(a+a \coth (e+f x))^3} \, dx}{6 f}+\frac{d \int \frac{1}{(a+a \coth (e+f x))^2} \, dx}{8 a f}\\ &=-\frac{d x^2}{16 a^3}+\frac{x (c+d x)}{8 a^3}-\frac{d}{36 f^2 (a+a \coth (e+f x))^3}-\frac{c+d x}{6 f (a+a \coth (e+f x))^3}-\frac{d}{32 a f^2 (a+a \coth (e+f x))^2}-\frac{c+d x}{8 a f (a+a \coth (e+f x))^2}-\frac{d}{16 f^2 \left (a^3+a^3 \coth (e+f x)\right )}-\frac{c+d x}{8 f \left (a^3+a^3 \coth (e+f x)\right )}+\frac{d \int 1 \, dx}{16 a^3 f}+\frac{d \int \frac{1}{a+a \coth (e+f x)} \, dx}{16 a^2 f}+\frac{d \int \frac{1}{(a+a \coth (e+f x))^2} \, dx}{12 a f}\\ &=\frac{d x}{16 a^3 f}-\frac{d x^2}{16 a^3}+\frac{x (c+d x)}{8 a^3}-\frac{d}{36 f^2 (a+a \coth (e+f x))^3}-\frac{c+d x}{6 f (a+a \coth (e+f x))^3}-\frac{5 d}{96 a f^2 (a+a \coth (e+f x))^2}-\frac{c+d x}{8 a f (a+a \coth (e+f x))^2}-\frac{3 d}{32 f^2 \left (a^3+a^3 \coth (e+f x)\right )}-\frac{c+d x}{8 f \left (a^3+a^3 \coth (e+f x)\right )}+\frac{d \int 1 \, dx}{32 a^3 f}+\frac{d \int \frac{1}{a+a \coth (e+f x)} \, dx}{24 a^2 f}\\ &=\frac{3 d x}{32 a^3 f}-\frac{d x^2}{16 a^3}+\frac{x (c+d x)}{8 a^3}-\frac{d}{36 f^2 (a+a \coth (e+f x))^3}-\frac{c+d x}{6 f (a+a \coth (e+f x))^3}-\frac{5 d}{96 a f^2 (a+a \coth (e+f x))^2}-\frac{c+d x}{8 a f (a+a \coth (e+f x))^2}-\frac{11 d}{96 f^2 \left (a^3+a^3 \coth (e+f x)\right )}-\frac{c+d x}{8 f \left (a^3+a^3 \coth (e+f x)\right )}+\frac{d \int 1 \, dx}{48 a^3 f}\\ &=\frac{11 d x}{96 a^3 f}-\frac{d x^2}{16 a^3}+\frac{x (c+d x)}{8 a^3}-\frac{d}{36 f^2 (a+a \coth (e+f x))^3}-\frac{c+d x}{6 f (a+a \coth (e+f x))^3}-\frac{5 d}{96 a f^2 (a+a \coth (e+f x))^2}-\frac{c+d x}{8 a f (a+a \coth (e+f x))^2}-\frac{11 d}{96 f^2 \left (a^3+a^3 \coth (e+f x)\right )}-\frac{c+d x}{8 f \left (a^3+a^3 \coth (e+f x)\right )}\\ \end{align*}

Mathematica [A]  time = 0.701283, size = 185, normalized size = 1.01 \[ \frac{\text{csch}^3(e+f x) \left (4 \left (6 c f (6 f x+1)+d \left (18 f^2 x^2+6 f x+1\right )\right ) \cosh (3 (e+f x))+27 (4 c f+d (4 f x+3)) \cosh (e+f x)+144 c f^2 x \sinh (3 (e+f x))+324 c f \sinh (e+f x)-24 c f \sinh (3 (e+f x))+72 d f^2 x^2 \sinh (3 (e+f x))+324 d f x \sinh (e+f x)-24 d f x \sinh (3 (e+f x))+135 d \sinh (e+f x)-4 d \sinh (3 (e+f x))\right )}{1152 a^3 f^2 (\coth (e+f x)+1)^3} \]

Antiderivative was successfully verified.

[In]

Integrate[(c + d*x)/(a + a*Coth[e + f*x])^3,x]

[Out]

(Csch[e + f*x]^3*(27*(4*c*f + d*(3 + 4*f*x))*Cosh[e + f*x] + 4*(6*c*f*(1 + 6*f*x) + d*(1 + 6*f*x + 18*f^2*x^2)
)*Cosh[3*(e + f*x)] + 135*d*Sinh[e + f*x] + 324*c*f*Sinh[e + f*x] + 324*d*f*x*Sinh[e + f*x] - 4*d*Sinh[3*(e +
f*x)] - 24*c*f*Sinh[3*(e + f*x)] - 24*d*f*x*Sinh[3*(e + f*x)] + 144*c*f^2*x*Sinh[3*(e + f*x)] + 72*d*f^2*x^2*S
inh[3*(e + f*x)]))/(1152*a^3*f^2*(1 + Coth[e + f*x])^3)

________________________________________________________________________________________

Maple [B]  time = 0.107, size = 791, normalized size = 4.3 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((d*x+c)/(a+a*coth(f*x+e))^3,x)

[Out]

1/a^3/f*(-4/f*d*(1/6*(f*x+e)*sinh(f*x+e)^3*cosh(f*x+e)^3-1/8*(f*x+e)*sinh(f*x+e)*cosh(f*x+e)^3+1/16*(f*x+e)*co
sh(f*x+e)*sinh(f*x+e)+1/32*(f*x+e)^2-1/36*sinh(f*x+e)^2*cosh(f*x+e)^4+13/288*sinh(f*x+e)^2*cosh(f*x+e)^2+1/72*
cosh(f*x+e)^2)+4/f*d*(1/6*(f*x+e)*sinh(f*x+e)^2*cosh(f*x+e)^4-1/12*(f*x+e)*sinh(f*x+e)^2*cosh(f*x+e)^2-1/12*(f
*x+e)*cosh(f*x+e)^2-1/36*sinh(f*x+e)*cosh(f*x+e)^5+1/36*cosh(f*x+e)^3*sinh(f*x+e)+1/24*cosh(f*x+e)*sinh(f*x+e)
+1/24*f*x+1/24*e)+1/f*d*(1/4*(f*x+e)*cosh(f*x+e)*sinh(f*x+e)^3-3/8*(f*x+e)*cosh(f*x+e)*sinh(f*x+e)+3/16*(f*x+e
)^2-1/16*sinh(f*x+e)^2*cosh(f*x+e)^2+1/4*cosh(f*x+e)^2)-3/f*d*(1/4*(f*x+e)*sinh(f*x+e)^2*cosh(f*x+e)^2-1/4*(f*
x+e)*cosh(f*x+e)^2-1/16*cosh(f*x+e)^3*sinh(f*x+e)+5/32*cosh(f*x+e)*sinh(f*x+e)+5/32*f*x+5/32*e)+4*d*e/f*(1/6*s
inh(f*x+e)^3*cosh(f*x+e)^3-1/8*cosh(f*x+e)^3*sinh(f*x+e)+1/16*cosh(f*x+e)*sinh(f*x+e)+1/16*f*x+1/16*e)-4*d*e/f
*(1/6*sinh(f*x+e)^2*cosh(f*x+e)^4-1/12*sinh(f*x+e)^2*cosh(f*x+e)^2-1/12*cosh(f*x+e)^2)-d*e/f*((1/4*sinh(f*x+e)
^3-3/8*sinh(f*x+e))*cosh(f*x+e)+3/8*f*x+3/8*e)+3*d*e/f*(1/4*sinh(f*x+e)^2*cosh(f*x+e)^2-1/4*cosh(f*x+e)^2)-4*c
*(1/6*sinh(f*x+e)^3*cosh(f*x+e)^3-1/8*cosh(f*x+e)^3*sinh(f*x+e)+1/16*cosh(f*x+e)*sinh(f*x+e)+1/16*f*x+1/16*e)+
4*c*(1/6*sinh(f*x+e)^2*cosh(f*x+e)^4-1/12*sinh(f*x+e)^2*cosh(f*x+e)^2-1/12*cosh(f*x+e)^2)+c*((1/4*sinh(f*x+e)^
3-3/8*sinh(f*x+e))*cosh(f*x+e)+3/8*f*x+3/8*e)-3*c*(1/4*sinh(f*x+e)^2*cosh(f*x+e)^2-1/4*cosh(f*x+e)^2))

________________________________________________________________________________________

Maxima [A]  time = 3.25172, size = 186, normalized size = 1.02 \begin{align*} \frac{1}{96} \, c{\left (\frac{12 \,{\left (f x + e\right )}}{a^{3} f} + \frac{18 \, e^{\left (-2 \, f x - 2 \, e\right )} - 9 \, e^{\left (-4 \, f x - 4 \, e\right )} + 2 \, e^{\left (-6 \, f x - 6 \, e\right )}}{a^{3} f}\right )} + \frac{{\left (72 \, f^{2} x^{2} e^{\left (6 \, e\right )} + 108 \,{\left (2 \, f x e^{\left (4 \, e\right )} + e^{\left (4 \, e\right )}\right )} e^{\left (-2 \, f x\right )} - 27 \,{\left (4 \, f x e^{\left (2 \, e\right )} + e^{\left (2 \, e\right )}\right )} e^{\left (-4 \, f x\right )} + 4 \,{\left (6 \, f x + 1\right )} e^{\left (-6 \, f x\right )}\right )} d e^{\left (-6 \, e\right )}}{1152 \, a^{3} f^{2}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((d*x+c)/(a+a*coth(f*x+e))^3,x, algorithm="maxima")

[Out]

1/96*c*(12*(f*x + e)/(a^3*f) + (18*e^(-2*f*x - 2*e) - 9*e^(-4*f*x - 4*e) + 2*e^(-6*f*x - 6*e))/(a^3*f)) + 1/11
52*(72*f^2*x^2*e^(6*e) + 108*(2*f*x*e^(4*e) + e^(4*e))*e^(-2*f*x) - 27*(4*f*x*e^(2*e) + e^(2*e))*e^(-4*f*x) +
4*(6*f*x + 1)*e^(-6*f*x))*d*e^(-6*e)/(a^3*f^2)

________________________________________________________________________________________

Fricas [A]  time = 2.09015, size = 709, normalized size = 3.87 \begin{align*} \frac{4 \,{\left (18 \, d f^{2} x^{2} + 6 \, c f + 6 \,{\left (6 \, c f^{2} + d f\right )} x + d\right )} \cosh \left (f x + e\right )^{3} + 12 \,{\left (18 \, d f^{2} x^{2} + 6 \, c f + 6 \,{\left (6 \, c f^{2} + d f\right )} x + d\right )} \cosh \left (f x + e\right ) \sinh \left (f x + e\right )^{2} + 4 \,{\left (18 \, d f^{2} x^{2} - 6 \, c f + 6 \,{\left (6 \, c f^{2} - d f\right )} x - d\right )} \sinh \left (f x + e\right )^{3} + 27 \,{\left (4 \, d f x + 4 \, c f + 3 \, d\right )} \cosh \left (f x + e\right ) + 3 \,{\left (108 \, d f x + 4 \,{\left (18 \, d f^{2} x^{2} - 6 \, c f + 6 \,{\left (6 \, c f^{2} - d f\right )} x - d\right )} \cosh \left (f x + e\right )^{2} + 108 \, c f + 45 \, d\right )} \sinh \left (f x + e\right )}{1152 \,{\left (a^{3} f^{2} \cosh \left (f x + e\right )^{3} + 3 \, a^{3} f^{2} \cosh \left (f x + e\right )^{2} \sinh \left (f x + e\right ) + 3 \, a^{3} f^{2} \cosh \left (f x + e\right ) \sinh \left (f x + e\right )^{2} + a^{3} f^{2} \sinh \left (f x + e\right )^{3}\right )}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((d*x+c)/(a+a*coth(f*x+e))^3,x, algorithm="fricas")

[Out]

1/1152*(4*(18*d*f^2*x^2 + 6*c*f + 6*(6*c*f^2 + d*f)*x + d)*cosh(f*x + e)^3 + 12*(18*d*f^2*x^2 + 6*c*f + 6*(6*c
*f^2 + d*f)*x + d)*cosh(f*x + e)*sinh(f*x + e)^2 + 4*(18*d*f^2*x^2 - 6*c*f + 6*(6*c*f^2 - d*f)*x - d)*sinh(f*x
 + e)^3 + 27*(4*d*f*x + 4*c*f + 3*d)*cosh(f*x + e) + 3*(108*d*f*x + 4*(18*d*f^2*x^2 - 6*c*f + 6*(6*c*f^2 - d*f
)*x - d)*cosh(f*x + e)^2 + 108*c*f + 45*d)*sinh(f*x + e))/(a^3*f^2*cosh(f*x + e)^3 + 3*a^3*f^2*cosh(f*x + e)^2
*sinh(f*x + e) + 3*a^3*f^2*cosh(f*x + e)*sinh(f*x + e)^2 + a^3*f^2*sinh(f*x + e)^3)

________________________________________________________________________________________

Sympy [A]  time = 4.70043, size = 1290, normalized size = 7.05 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((d*x+c)/(a+a*coth(f*x+e))**3,x)

[Out]

Piecewise((36*c*f**2*x*tanh(e + f*x)**3/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864
*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 108*c*f**2*x*tanh(e + f*x)**2/(288*a**3*f**2*tanh(e + f*x)**3 + 86
4*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 108*c*f**2*x*tanh(e + f*x)/(288*
a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 3
6*c*f**2*x/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 28
8*a**3*f**2) - 108*c*f*tanh(e + f*x)**3/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864
*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) - 72*c*f*tanh(e + f*x)**2/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3
*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 12*c*f/(288*a**3*f**2*tanh(e + f*x)**3
 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 18*d*f**2*x**2*tanh(e + f*x
)**3/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3
*f**2) + 54*d*f**2*x**2*tanh(e + f*x)**2/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 86
4*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 54*d*f**2*x**2*tanh(e + f*x)/(288*a**3*f**2*tanh(e + f*x)**3 + 86
4*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 18*d*f**2*x**2/(288*a**3*f**2*ta
nh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) - 87*d*f*x*tanh
(e + f*x)**3/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) +
288*a**3*f**2) - 9*d*f*x*tanh(e + f*x)**2/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 8
64*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 63*d*f*x*tanh(e + f*x)/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**
3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) + 33*d*f*x/(288*a**3*f**2*tanh(e + f*x)
**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2) - 45*d*tanh(e + f*x)**3/(2
88*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(e + f*x) + 288*a**3*f**2)
- 48*d*tanh(e + f*x)**2/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f**2*tanh(
e + f*x) + 288*a**3*f**2) + 11*d/(288*a**3*f**2*tanh(e + f*x)**3 + 864*a**3*f**2*tanh(e + f*x)**2 + 864*a**3*f
**2*tanh(e + f*x) + 288*a**3*f**2), Ne(f, 0)), ((c*x + d*x**2/2)/(a*coth(e) + a)**3, True))

________________________________________________________________________________________

Giac [A]  time = 1.16859, size = 204, normalized size = 1.11 \begin{align*} \frac{{\left (72 \, d f^{2} x^{2} e^{\left (6 \, f x + 6 \, e\right )} + 144 \, c f^{2} x e^{\left (6 \, f x + 6 \, e\right )} + 216 \, d f x e^{\left (4 \, f x + 4 \, e\right )} - 108 \, d f x e^{\left (2 \, f x + 2 \, e\right )} + 24 \, d f x + 216 \, c f e^{\left (4 \, f x + 4 \, e\right )} - 108 \, c f e^{\left (2 \, f x + 2 \, e\right )} + 24 \, c f + 108 \, d e^{\left (4 \, f x + 4 \, e\right )} - 27 \, d e^{\left (2 \, f x + 2 \, e\right )} + 4 \, d\right )} e^{\left (-6 \, f x - 6 \, e\right )}}{1152 \, a^{3} f^{2}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((d*x+c)/(a+a*coth(f*x+e))^3,x, algorithm="giac")

[Out]

1/1152*(72*d*f^2*x^2*e^(6*f*x + 6*e) + 144*c*f^2*x*e^(6*f*x + 6*e) + 216*d*f*x*e^(4*f*x + 4*e) - 108*d*f*x*e^(
2*f*x + 2*e) + 24*d*f*x + 216*c*f*e^(4*f*x + 4*e) - 108*c*f*e^(2*f*x + 2*e) + 24*c*f + 108*d*e^(4*f*x + 4*e) -
 27*d*e^(2*f*x + 2*e) + 4*d)*e^(-6*f*x - 6*e)/(a^3*f^2)